arXiv:1501.01753vl [cond-mat.mtrl-sci] 8 Jan 2015 


Using ultrashort optical pulses to couple ferroelectric and ferromagnetic order in an 

oxide heterostructure 

Y. M. Sheu, 1 S. A. Trugman, 1 L. Yan, 1 Q. X. Jia, 1 A. J. Taylor, 1 and R. P. Prasankumar 1 

1 Center for Integrated Nanotechnologies, MS K771, 

Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA 
(Dated: January 9, 2015) 

A new approach to all-optical detection and control of the coupling between electric and magnetic 
order on ultrafast timescales is achieved using time-resolved second harmonic generation (SHG) to 
study a ferroelectric (FE)/ferromagnet (FM) oxide heterostructure. We use femtosecond optical 
pulses to modify the spin alignment in a Bao.iSro. 9 Ti 03 (BSTO)/Lao. 7 Cao. 3 MnC >3 (LCMO) het¬ 
erostructure and selectively probe the ferroelectric response using SHG. In this heterostructure, 
the pump pulses photoexcite non-equilibrium quasiparticles in LCMO, which rapidly interact with 
phonons before undergoing spin-lattice relaxation on a timescale of tens of picoseconds. This re¬ 
duces the spin-spin interactions in LCMO, applying stress on BSTO through magnetostriction. This 
then modifies the FE polarization through the piezoelectric effect, on a timescale much faster than 
laser-induced heat diffusion from LCMO to BSTO. We have thus demonstrated an ultrafast indi¬ 
rect magnetoelectric effect in a FE/FM heterostructure mediated through elastic coupling, with a 
timescale primarily governed by spin-lattice relaxation in the FM layer. 


Magnetoelectric multiferroics have attracted much re¬ 
cent attention, largely due to the scarcity of naturally 
occurring materials in which magnetic and electric or¬ 
der coexist [1-4]. However, the microscopic mechanisms 
underlying magnetoelectric (ME) coupling between these 
order parameters have been difficult to unravel, hinder¬ 
ing the development of single phase multiferroics with 
strong ME coupling at useful temperatures [1, 2, 5]. Ar¬ 
tificial multiferroic composites [6-9] offer an appealing 
alternative in which the ME coupling can be engineered 
through proper choice of the interface geometry and con¬ 
stituent materials [4]. The most well known approach for 
achieving strong ME coupling in these systems uses strain 
to indirectly couple ferroelectric (FE) and ferromagnetic 
(FM) order in two-phase layered composites [5, 7, 10]. 
In these composites, the ME effect is essentially a prod¬ 
uct of the magnetostrictive effect in the FM layer and the 
piezoelectric effect in the FE layer [5, 10-13]. A magnetic 
(B) held is used to modify the spin alignment in the FM 
layer, which changes the lattice constant through mag¬ 
neto striction. This in turn applies stress on the FE layer 
and modifies the FE polarization through the piezoelec¬ 
tric effect [12-14], significantly increasing the ME cou¬ 
pling as compared to single-phase multiferroics [5, 10]. 
Several reviews describing the ME effect in these com¬ 
posites have been recently published [5, 10, 12, 15]. 

Despite these impressive advances, an important as¬ 
pect of multiferroics has received relatively little atten¬ 
tion: namely, their dynamic properties [2, 3, 10], which 
will impact many of their potential applications. In fact, 
the ultimate timescales limiting the speed of ME coupling 
in these materials, as well as the physical processses gov¬ 
erning them, remain essentially unexplored. Femtosec¬ 
ond optical pulses are particularly useful in this regard 
[16-31], since they can be used to both interrogate and 
control the ME response in a non-contact manner on ul¬ 


trafast timescales, much faster than by directly applying 
magnetic or electric fields. Furthermore, specific order 
parameters can be directly accessed using time-resolved 
second harmonic generation (TR-SHG), making it espe¬ 
cially attractive for studying both magnetically [20] and 
orbitally ordered [32] materials as well as multiferroics 

[18] and interfacial properties of novel layered materials 

[19] in the time domain. 

Here, we use TR-SHG to explore and optically 
manipulate the coupling between ferroelectric and 
ferromagnetic order in an oxide heterostructure for the 
first time, inspired by the strain/stress-based approach 
described above. By separating the different contri¬ 
butions to the TR-SHG signal in the time domain, 
we discovered that the timescale dominating the ME 
response is governed by demagnetization of the FM 
layer through spin-lattice relaxation. More specifi¬ 
cally, optically perturbing magnetic order in the FM 
layer imposes lateral stress on the FE layer through 
magnetostriction, modifying FE order within tens of 
picoseconds (ps). This demonstrates that femtosecond 
optical pulses can be used not only to give insight into 
the microscopic mechanisms underlying ME coupling in 
complex oxide heterostructures, but also to manipulate 
the ME response in these systems on ultrafast timescales. 

Results 

Static and time-resolved second harmonic gen¬ 
eration (SHG) experiments. The heterostructure 
studied here consists of a 50 nm thick film of ferroelectric 
BSTO and a 50 nm thick film of ferromagnetic LCMO, 
grown on a MgO substrate; more detail on sample 
fabrication is given in Methods and in Ref. [33]. We 
chose BSTO since lattice strain introduced through 
epitaxial growth can enhance ferroelectricity and the FE 
transition temperature (Tfe), as theoretically predicted 
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FIG. 1. Experimental setup, static SHG characteri¬ 
zation, and piezoelectric response to in-plane stress. 

(a) Experimental schematic for TR-SHG in reflection. 9 is the 
angle between the sample normal and the light propagation 
direction. <ji is the incident polarization of the fundamental 
light, which is controlled by a half-wave (A/2) plate. Inci¬ 
dent p and s polarizations correspond to (j>= 0° (180°) and 
90° (270°) (b) Polar plot of the measured SHG intensity at 
10 K. The blue and red symbols are the p and s polarized 
SHG signals, I p and I s , respectively, plotted as a function of 
the incident light polarization. The cyan dashed and green 
solid lines are numerical fits to I v and I s with d 31 =-l, d 15 =- 
1.2, and d 3 3=7. (c) Schematic diagram of piezoelectric re¬ 
sponse along z to applied in-plane stress {e xx and e yy ). (d) 
Static temperature-dependent SHG for (p,p) and (s,p) polar¬ 
ization combinations. In LCMO, Tq ~240 K, while in BSTO, 
Tfe ~215 K [33], as indicated by the arrow. 


[34, 35] and experimentally observed [33]. This coupling 
between strain and FE order makes it a good candidate 
for inducing indirect ME effects through elastic coupling 
with a magnetic material. Similarly, LCMO was chosen 
since a strong ME effect was previously obtained through 
magnetostriction under an applied DC magnetic field 
when it was incorporated into a bilayer structure with 
lead zirconate titanate (PZT) [5, 14, 36]. The concept 
behind our experiment is thus to demagnetize LCMO 
with femtosecond optical pulses, modifying its lattice 
constant through magnetostriction; the resulting stress 
on BSTO activates its piezoelectric response, changing 
FE order and indirectly inducing ME coupling within 
the demagnetization time. 

Photoinduced changes in the static FE polarization of 
BSTO, P s , can be measured using SHG (Figs, la, lb, and 
Id), which is well known to directly probe FE order [37- 
40]. The SHG polarization generated by the fundamen¬ 
tal light fields Ej and £/ is Pi(2uj) = eodijkEj(oj)Ek(w), 
where w is the frequency and dijk is the nonlinear opti¬ 
cal susceptibility tensor [39-41], Our BSTO films have 
Giv symmetry and a FE polarization along the sample 
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FIG. 2. Photoinduced coupling between ferroelectric 
and ferromagnetic order, (a) Schematic diagram of the 
interplay between the lattice, the magnetization M, and the 
FE polarization P at various pump-probe delays. After pho¬ 
toexcitation, thermal stress on BSTO (thin red arrows) ini¬ 
tially results from e — ph coupling in LCMO within 1 ps. Its 
magnitude gradually increases after relaxation of the lattice 
contraction in LCMO through magnetostriction (thin blue ar¬ 
rows) on a timescale of ~ 50 ps, causing P to decrease. Ther¬ 
mal diffusion from LCMO to both BSTO and the substrate 
then takes place at longer timescales, t >0.5 ns, finally dissi¬ 
pating through the substrate, (b) Time-resolved p-polarized 
SHG as a function of the incident light polarization at 10 K. 
P 0 -P 3 are the FE polarization responses at different times to 
changes in the magnetization (arising from M 0 -M 3 in (a)). 


normal (z in Fig. la) due to the in-plane compressive 
strain induced by lattice matching to LCMO [33]. There 
are only three non-vanishing tensor components for this 
crystal symmetry: di 5 ,<i 3 i, and d 33 . By varying the po¬ 
larization of the incident fundamental and detected SHG 
light (Fig. la) while measuring the SHG intensity I (Fig. 
lb), these different components can be extracted, as dis¬ 
cussed in Methods. 

To optically explore the ME response of our 
BSTO/LCMO heterostructure, we use the majority 
of the output laser intensity to photoexcite LCMO (at 
1.59 eV, far below the BSTO band gap of ^3.3-3.6 eV) 
[33]) and the remainder to measure the resulting changes 
in P s using SHG. Photoexcited quasiparticles (unbound 
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FIG. 3. Evidence for a non-equilibrium, nonlinear ori¬ 
gin of the TR-SHG response, (a) Polarization-dependent 
changes in I p measured at 10 K for various time delays. The 
data is plotted as symbols and numerical calculations for t=- 
10, 50 ps are plotted as dashed lines. Before t = 0, the 
change is due to long-lived residual heating from LCMO to 
BSTO. After t = 0 the change arises from reductions in cfai 
and dss (more detail is given in Methods), (b) Polarization- 
dependent changes in I p measured at t=200 ps for different 
temperatures. 


electrons and holes) in LCMO will relax through 
different processes that have been extensively studied 
using pump-probe spectroscopy [42-44]. Initially, 
the non-equilibrium quasiparticles rapidly lose energy 
through electron-phonon (e — ph) coupling (within ^1 
ps), after which the lattice exchanges energy with the 
spins through spin-lattice (s — l) coupling on a timescale 
of tens to hundreds of ps. On longer timescales (>1 
nanosecond (ns)) LCMO returns to equilibrium as the 
remaining energy is dissipated through the substrate. 
More detail on quasiparticle dynamics in manganites is 
provided in Supplementary Note 1. 

Data analysis. We can relate this to the timescales 
on which the FE polarization changes in our TR-SHG 
experiments (Fig. 2). Fig. 2a schematically shows the 


different processes occurring for different pump-probe de¬ 
lays t. Before t = 0, a small negative change in the 
p-polarized SHG intensity (A I p ) is observed in all polar 
combinations at low temperatures (Pq), corresponding to 
residual heating that is not completely dissipated before 
the arrival of the next pump pulse. This change in I p , 
where the SHG intensity decreases for all incident light 
polarizations, as more clearly seen in Figs. 2b and 3a 
for negative time delays, is equivalent to increasing the 
sample temperature in the SHG polar plot of Fig. lb (or 
the temperature-dependent plots in Fig. Id), since be¬ 
fore t =0 the electron, spin, and lattice subsystems share 
a common temperature. 

In contrast, at early times (t ^1-2 ps), we observe a 
polarization-dependent change in I p (Pi), where the SHG 
intensity both increases and decreases for different polar¬ 
izations, as shown in Figs. 2b and Fig. 3. This cannot be 
obtained by simply changing the sample temperature; it 
is a non-equilibrium effect, occurring before the different 
subsystems can reach a common temperature. This is 
due to the elevated temperature of the lattice in LCMO 
after e — ph coupling, causing it to expand, which in turn 
initiates an in-plane thermal stress e xx and £ yy on BSTO 
(without changing its temperature) (Fig. 2a). This then 
initiates an out-of-plane piezoelectric response through 
Dt = eijkSjk , where Di is the dielectric displacement 
and eijk is the piezoelectric coefficient. In the Chv sym¬ 
metry, e 3 i is non-zero (more detail is given in Methods), 
linking e xx and e yy to D z (Fig. lc) and leading to a 
photoinduced change in the SHG intensity. It takes ^7 
ps for the strain induced by in-plane stress to propagate 
through the BSTO layer (estimated from the speed of 
sound in BSTO), which agrees relatively well with the 
timescale of the initial rise in A/ p , as shown in Fig. 4 
(i.e. there is no rapid ~1 ps rise time, unlike the pho¬ 
toinduced change in reflectivity, A P/P). However, the 
maximum change in I p (represented by P 2 ) does not oc¬ 
cur on this timescale, but on a much slower timescale of 
~50-100 ps (Figs. 2b, and 3a). 

We can gain insight into the origin of this large change 
in I p (and thus P s ) by first noting that it primarily 
occurs below the FM ordering temperature, Tq ~240 
K (Fig. 3b), indicating that FM order in LCMO is 
linked to the observed A I p . In addition, the timescale 
is comparable to that for s — l relaxation in LCMO, 
as determined through comparison of our time-resolved 
SHG data to optical pump-probe measurements on our 
BSTO/LCMO heterostructure (Fig. 4). This suggests 
that the observed photoinduced ME response is initiated 
by a change in the lattice constant through magnetostric¬ 
tion [45], which arises from photoinduced demagnetiza¬ 
tion of LCMO. More specifically, s — l relaxation reduces 
the magnetization in the FM ordered state, gradually 
decreasing the magnitude of the existing in-plane lattice 
contraction due to spin-spin interactions (since the spins 
in LCMO are aligned in-plane) and applying stress to 
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FIG. 4. Comparison of A R/R and A I/I at (a) 10 
K and (b) above 250 K. AR/R is measured with a p- 
polarized probe and AI/I is taken in the ( p,p ) configuration. 
The solid lines are the experimental data and the dashed lines 
represent a double exponential fit. Additional discussion of 
the comparison between the A R/R and AI/I data is given 
in Methods and Supplementary Note 1. 


BSTO (Fig. 2a, t ~ 50 ps); more detail on general mag¬ 
netostriction is given in Supplementary Note 2. Then, as 
for the early time dynamics, the resulting increase in the 
in-plane tensile stress on BSTO initiates a piezoelectric 
effect along the z direction (Fig. lc), reducing the dielec¬ 
tric displacement D z and modifying the FE polarization 
P s (Fig. 2a). 

Additional experimental evidence enables us to rule 
out other possibilities that could lead to the observed 
change in A I p . First, A I p at low temperatures does not 
originate from a photoinduced modulation in the reflec¬ 
tivity due to changes in the optical constants; A R/R is 
an order of magnitude smaller than A I // (the normal¬ 
ized photoinduced change in SHG intensity for a given 
polarization configuration) at low temperatures (Fig. 4), 
and unlike A I // , it does not vary with the incident probe 
polarization (Supplementary Figure 1). Second, there is 
no observable SHG signal, either static ( I p ) or photoin¬ 
duced (AI/I), from a single LCMO film, verifying that 


A I p in our BSTO/LCMO heterostructure does not re¬ 
sult from SHG in LCMO. Third, there is no measurable 
AI/I signal from a single BSTO film upon 1.59 eV exci¬ 
tation, since the pump is far below the gap [33]. Fourth, 
we note that the observed features are independent of the 
pump polarization, consistent with photoinduced demag¬ 
netization in LCMO. Fifth, Kerr rotation of the reflected 
light at both 1.55 eV and 3.1 eV, induced by demagneti¬ 
zation in LCMO, does not contribute to A I p . Previous 
work has shown that the Kerr rotation in manganites for 
similar pump fluences is ~0.01-0.1 degrees [44, 46]. This 
small rotation is unlikely to result in the large changes in 
SHG intensity that we observed. In addition, if it were 
to contribute significantly to AI/I , it would also be ob¬ 
served in TR-SHG measurements on single LCMO films, 
where Kerr rotation dominates polarization-dependent 
measurements [44, 46]. Moreover, in the unlikely case 
that Kerr rotation in this heterostructure is larger than 
that of a single LCMO layer upon photoexcitation, the 
rotated fundamental or SHG polarizations would induce 
phase shifts in our polarization-dependent TR-SHG mea¬ 
surements, which are not observed (Fig. 2b and Fig. 3). 

Finally, any contributions from surface SHG, including 
any potential time-dependent interference with the con¬ 
tribution from FE order, are expected to be negligible. 
The temperature-dependent static SHG signal (Fig. Id) 
originates primarily from bulk FE order in BSTO [33]. 
Surface and interface contributions from this layer have 
the same symmetry, C 4 V , as the bulk FE order. Any 
interference between these contributions will be either 
in-phase or 180° out-of-phase, but will remain constant 
over the entire sample, as a spatially homogeneous SHG 
intensity is observed after SHG saturation [33]. Well 
above Tfe, the residual SHG signal likely originates from 
the surface or interface and mixes with other potential 
contributions, such as slow carrier transfer [33] and size 
effects [37, 47] in thin films. Upon photoexcitation, 
surface SHG does not contribute to the AI/I signal, 
since the BSTO surface is not excited. Furthermore, we 
can exclude any contribution from changes in interfacial 
SHG due to demagnetization of LCMO, because no 
TR-SHG signal was obtained from the interface between 
LCMO and the MgO substrate. These considerations 
allow us to conclude that the AI/I signal from our 
heterostructure originates from changes in P s that are 
initiated by photoinduced changes in LCMO. 

Discussion 

It is worth re-emphasizing that, as discussed in the intro¬ 
duction, the photoinduced ME response measured here 
is analogous to the ME effect observed when an external 
B field is applied to enhance the lattice contraction 
through magnetostriction in a FE/FM bilayer, inducing 
a piezoelectric response in the ferroelectric layer [14]. 
The major advantage of our approach lies in the tempo¬ 
ral degree of freedom offered by femtosecond laser pulses, 
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enabling us to optically manipulate the ME response in 
the time domain and also discriminate between different 
contributions to the ME response. For example, we can 
temporally separate the comparatively small changes in 
P s at early times due to heating-induced expansion of 
LCMO (Fig. 2, t ~1 ps) from the larger changes in P s 
on longer timescales due to the reduction in the LCMO 
lattice contraction through magnetostriction (Fig. 2, 
t ~50 ps). The difference in the magnitude of these 
changes is likely due to the strong interaction between 
the ordered spins in LCMO (causing the lattice to 
contract through magnetostriction), which only permits 
a small amount of heating-induced thermal stress at 
early times. As a result, the in-plane tensile stress does 
not significantly increase before s — l relaxation takes 
place. Therefore, strain relaxation throughout the entire 
BSTO film only occurs after tens of ps, setting the 
timescale for the maximum changes in FE polarization 
(Fig. 3a). Finally, we note that at long (hundreds of 
ps to ns) timescales, the remaining changes in the FE 
polarization are due to thermal diffusion from LCMO to 
BSTO. 

In conclusion, we used femtosecond optical pulses 
to manipulate the coupling between FE and FM order 
in a complex oxide heterostructure, employing time- 
resolved SHG to optically perturb the magnetization 
in LCMO and selectively probe the associated change 
in the ferroelectric properties of BSTO. We find that 
although lattice heating in LCMO occurs within ~1 ps, 
it does not significantly apply lateral stress on BSTO 
because of the strong lattice contraction originating from 
spin-spin interactions in LCMO. This lattice contraction 
is reduced through spin-lattice relaxation, increasing 
the tensile stress on BSTO through magnetostriction 
in LCMO and leading to a piezoelectric response that 
changes the FE polarization. Therefore, coupling 
between FE and FM order in BSTO/LCMO is induced 
within tens of picoseconds, mediated through elastic 
coupling between the BSTO and LCMO layers. Fem¬ 
tosecond optical spectroscopy can thus shed light on the 
mechanisms underlying ME coupling in complex oxide 
heterostructures, while presenting intriguing possibilities 
for future high speed optically controlled magnetoelec¬ 
tric devices (e.g., through coherent spin control [48-50], 
which would enable faster manipulation of spin dynamics 
and minimize undesired heat flow in device applications.) 

Methods 

Sample fabrication and characterization. The 

samples used in our studies are 50 nm thick BSTO 
films grown by pulsed laser deposition on 50 nm thick 
optimally doped manganite films (Lao. 7 Cao. 3 Mn 03 ), 
using (100) MgO substrates. The substrate temperature 
during film growth is initially optimized and maintained 
at 750 °C. The oxygen pressure during deposition is 100 
mTorr. The samples are cooled to room temperature in 


pure oxygen (at a pressure of 350 Torr) without further 
thermal treatment. X-ray diffraction measurements 
reveal a compressive strain in epitaxially grown BSTO, 
while LCMO on MgO is relaxed. This causes an increase 
in tetragonality, even at room temperature, enhancing 
the FE Tq to ~215 K, which is significantly higher than 
the bulk Tq ~ 75 K. More detailed characterization of 
these samples is discussed in Ref. [33]. 

TR-SHG experimental details. Our time-resolved 
SHG experiments are based on an amplified Ti:sapphirc 
laser system producing pulses at a 250 kHz repetition 
rate with ~100 fs duration and energies of ~4 yj at a 
center wavelength of 780 nm (1.59 eV). SHG is generated 
from BSTO at 3.18 eV using a 1.59 eV probe beam, 
after which it is detected by a photomultiplier tube 
(PMT) using lock-in detection after filtering out the 
fundamental signal. The fundamental light polarization 
is controlled by a half-wave plate (A/2 in Fig. la) and 
the SHG signal is detected for either p or s polarizations. 
The standard laser fluence used for SHG in our mea¬ 
surements is Fq ^0.25 mJ cm” 2 and the pump fluence is 
~1 mJ cm -2 , creating a carrier density of ~10 2 ° cm -3 
(~T0 -2 per unit cell). In the C± v symmetry, SHG can 
be generated and probed when the fundamental light 
is incident at an angle away from the sample normal, 
thus exhibiting a sinusoidal dependence on the inci¬ 
dent angle 0 (as observed in separate angle-dependent 
measurements on our samples). Therefore, we used a 
reflection geometry in which the fundamental light was 
incident at a ~20° angle to obtain the data shown in 
our manuscript. The pump beam was incident at a 
smaller angle of ~10-15°, and its specular reflection was 
completely separated from the probe beam path. The 
temporal delay is achieved by a mechanical delay stage 
that allows us to vary the probe path length. Finally, 
ref. [33] gives a detailed static SHG characterization of 
these samples. Here, all time-resolved data was taken 
after saturation of the static SHG signal [33]. 

SHG analysis. For the third rank tensor dijk, 
the indices jk can be replaced by Z=l,2,...,6, corre¬ 
sponding to xx, yy, zz, yz, xz, and xy, using simple 
symmetry arguments that reduce 27 tensor components 
to 18. Further reduction can be achieved through noting 
that the C^ v symmetry has d 3 \ = d zxx = d 32 = d zyy , 
di 5 = d xxz = d 2 4 = d yyz and d 33 = d zzz . Therefore, 
we use the three independent indices 15, 31, and 33 
here. To represent the fourth rank stress tensor e^fe, 
the first i= 1, 2. ..,6 as defined above, and j,k= 1, 2, 3 
correspond to the Cartesian coordinate x,y,z. It also 
can be reduced according to the above discussion of 
indexes j, k. 

Fig. lb depicts the static SHG signals for p and s 
polarizations, I p and I s . In reflection (Fig. la), cf> de¬ 
termines the polarization of the fundamental light that 
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generates SHG. With a unit E field of the fundamen¬ 
tal light incident at an angle 6 (Fig. la), we have 
components (. E x {oj),E y (u)), E 2 (w))=(cos</>cos0,sin<(>,- 

cos^sinf?), where 0=0 is along 2 and <^=0° and 90° 
represent p and s polarizations of the incident funda¬ 
mental light, respectively. When the detection polarizer 
is set to p or s, we detect the SHG dipole components 
(P x (2u),P y (2u),P z (2u)) along (1,0,1) and (0,1,0), respec¬ 
tively, with each component linked to (di 5 = ^ 24,^31 = 

d 3 2 ,d 33 ) by 

P x ( 2w) = 2e 0 di 5 E x (uj)E z (u}), (1) 
P y (2u}) = 2e 0 di 5 E y (uj)E z (uj), (2) 
P z (2w) = eo(d 31 £g(w) + d 33 Ey(u)) + d 33 E 2 (oj)). (3) 

We thus can derive the SHG intensities I p and I s as 
functions of </> (Fig. lb): 

I p = (a cos 2 (j) + b sin 2 (j )) 2 , (4) 

I s = (csin(2^)) 2 , (5) 

where a, b and c are functions of 6 and linear combina¬ 
tions of di 5 , d 3 \, and ^ 33 . This enables us to detect dif¬ 
ferent components of dijk through different polarization 
combinations of (incident light, SHG). The combination 
of (45°,s) probes c, associated with d\ 3 , while (s,p) de¬ 
tects 6 , associated with d 3 i. The (p,p) configuration mea¬ 
sures a as a linear combination of all three independent 
components of d^k, allowing us to extract <^33 (Supple¬ 
mentary Table 1) from the measurements in the (45°,s) 
and (s,p) configurations. 

We can use the above analysis methods to numerically 
fit our TR-SHG data (Fig. 3a). To do this, we need 
to know the new values of a and b after photoexcita¬ 
tion, which can then be inserted into Eq. (4). First, 
we note that the polarization dependence of A I p for 
t > 0 (i.e., A I p decreases for (s,p) and increases for 
(p,p)) (Fig. 3) implies that a decrease in b is accom¬ 
panied by an increase in a. The TR-SHG signal in the 
(s,p) configuration is solely due to b (originating from g? 3 i 
(E q. (4))),which has the opposite sign of a (Supplemen¬ 
tary Table 1). Therefore, the photoinduced decrease in 
d 3 i implies an increase in a, as observed for the (p, p) con¬ 
figuration in Fig. 3. Next, we consider the coefficient c, 
which solely originates from ^ 15 . There are three signifi¬ 
cant features we notice from the photoinduced change in 
^ 15 . First, A I /1 is much smaller than that measured for 
the (p,p) and (s,p) configurations (Supplementary Figure 
1). Second, the change is comparable to the modulation 
arising from the change in optical constants (Supplemen¬ 
tary Figure lc and Fig. 4). Third, it does not have sig¬ 
nificant temperature dependence, unlike the other two 
polar combinations (Supplementary Figure 1). 

We conclude that d\ 3 has an insignificant contribution 
to a, and the photoinduced changes in d 15 could be 


strongly influenced by photoinduced reflectivity changes 
in BSTO. These considerations allowed us to simulate 
the photoinduced changes in d 33 and ^ 31 , revealing that 
Ae? 33 / d 33 ~ —4% and Ad 3 i/d 3 i ~ — 6 %. We thus can 
use these reductions in d 3 \ and d 33 , together with the 
isotropic reduction of d 3 \,d 33 and G?i 5 due to residual 
heating, to numerically fit A/p (dashed lines in Fig. 3a). 
These transient changes lead to the feature observed in 
I p after to=0 (Figs. 3a and Fig. 3b). 
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